Angular size and emission time scales of relativistic fireballs 
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ABSTRACT 

The detection of delayed X-ray, optical and radio emission, "afterglow," associated 
with 7-ray bursts (GRBs) is consistent with models, where the bursts are produced 
by relativistic expanding blast waves, driven by expanding fireballs at cosmological 
distances. In particular, the time scales over which radiation is observed at different 
wave bands agree with model predictions. It had recently been claimed that 
the commonly used relation between observation time t and blast wave radius 
r, t = r/2 r y 2 (r)c where 7 is the fluid Lorentz factor, should be replaced with 
t = r/167 2 (r)c due to blast wave deceleration. Applying the suggested deceleration 
modification would make it difficult to reconcile observed time scales with model 
predictions. It would also imply an apparent source size which is too large to allow 
attributing observed radio variability to diffractive scintillation. We present a detailed 
analysis of the implications of the relativistic hydrodynamics of expanding blast waves 
to the observed afterglow. We find that modifications due to shock deceleration are 
small, therefore allowing for both the observed afterglow time scales and for diffractive 
scintillation. We show that at time t the fireball appears on the sky as a narrow ring of 
radius h = r /^{r) and width Ah/h ~ 0.1, where r and t are related by t = r /2j 2 (r)c. 



Subject headings: gamma rays: bursts 



Introduction 



The availability of accurate positions for GRBs from the BeppoSAX satellite flCosta et al. 



19974 Feroci et al. 1997, 


Heise et al. 1997 


Costa et al. 1997c) allowed for the first time to detect 


delayed emission associated with GRBs in X-ray (posta et al. 1997a, b,Feroci et al. 1997 


Pirc 


|et al. 19974b), optical (Groot et al. 1997 


, Sahu et al. 1997], |van Paradijs et al. 1997 


. I 


3ond 


1997, Djorgovski et al. 1997) and radio (F 


rail & Kulkarni 1997) wave-bands. The detection of 



absorption lines in the optical afterglow of GRB970508 provided the first direct estimate of source 
distance, constraining the redshift of GRB970508 to 0.8 < z < 2.3 ( [Metzger et al. 1997] ) . Observed 
X-ray to radio afterglows are broadly consistent with models based on relativistic blast waves at 



cosmological distances ( Paczyhski fc Rhoads 1993 , Meszaros fc Rees 1997 , Vietri 1997 , Waxman 



1997a ,b, Wijers, Rees fc Meszaros 1997 ). Using these models, combined radio and optical data 
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allowed for the first time to directly estimate the total GRB energy, implying an energy of 
~ 10 52 erg for GRB970508 flWaxman 1997% 



In fireball models of GRB afterglow, the energy released by an explosion, ~ 10 52 erg, is 
converted to kinetic energy of a thin baryon shell expanding at ultra-relativistic speed. After 
producing the GRB, the shell impacts on surrounding gas, driving an ultra-relativistic shock into 
the ambient medium. After a short transition phase, the expanding blast wave approaches a 
self-similar behavior where the expansion Lorentz factor decreases with radius as 7 oc r~ 3 / 2 . The 
expanding shock continuously heats fresh gas and accelerates relativistic electrons, which produce 
the observed radiation through synchrotron emission. 

Photons emitted at radius r with frequency u, measured in the shell rest frame, are observed 
over a wide frequency range, 2^{r^)v to v/2 r y(r), and wide time range, since photons emitted further 
from the source-observer line of sight arrive at later time and are observed to have lower energy. 
There is no unique relation, therefore, between the radius r and the time or frequency at which 
radiation emitted at r is observed. However, since most of the emission from radius r detected 
by a distant observer originates from a disk of radius ~ r/~f{r) around the source-observer line of 
sight, and since photons emitted at a distance r/j(r) from the line of sight with frequency v (in 
the shell frame) are observed with frequency ^{r)v and are delayed (compared to photons emitted 
on the line of sight) by r/27 2 (r)c, it is commonly assumed that at observed time t = r/2^ 2 (r)c the 
flux peaks at frequency v' = 7(r)z^(r), where v(r) is the frequency at which the emission peaks 
in the shell rest frame. Based on the same argument, it is assumed that photons of frequency 
v 1 = 7(r)^(r) arrive over a time scale ~ r /2^ 2 (r)c. 



It had recently been claimed ( Sari 19971) that due to shock deceleration the commonly used 



relation t = r/2^ 2 (r)c should be replaced by t = r/l6~/ 2 (r)c. Applying the suggested deceleration 
modification would make it difficult to reconcile observed afterglow time scales with blast wave 
model predictions. This modification may have other implications. For example, shortly after 
the suggestion by Goodman (1997), that local inter-stellar scintillation may modulate the radio 
flux of GRB afterglows, radio variability on time scale consistent with scintillation origin has 
been observed ( [Frail et al. 1997 ). However, while the source size based on using the relation 



t = r/167 2 (r)c is ~ 4 times larger than required to allow for diffractive scintillation, implying 
that only refractive scintillation is likely (Goodman 1997), the observed large amplitude of radio 
variability suggests that modulation, if due to scintillation, should be produced by diffraction. 
Using the relation t = r /2^f 2 (r)c, on the other hand, implies a source size which is compatible 
with that required for diffractive scintillation. Similarly, the larger size implied by the suggested 
modification implies that it would be difficult to observe microlensing of the optical emission 
(|Loeb & Perna 1997] ). 



In §2 we present a detailed analysis of the implications of the relativistic hydrodynamics of 
expanding blast waves to the observed afterglow. We find that the modification due to shock 
deceleration of the numerical coefficient in the expression t = r/27 2 (r)c, giving the time at which 
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the flux peaks at observed frequency v' = j(r)u(r), is small. It should be emphasized, that the 
exact value of the numerical coefficient depends not only on the hydrodynamics, but also on 
the details of the electron energy distribution and on the magnetic field distribution behind the 
shock. Due to the lack of a theory describing these distributions, the exact value of the numerical 
coefficient can not be determined. However, it can be shown that the effect of deceleration is small. 
This is due to the fact that the relation t = r/\6j 2 (r)c gives the delay of photons emitted from 
the fireball at radius r from a point on the shock front on the line of sight, while most photons 
suffer longer delays, since they are emitted from a shell of finite thickness behind the shock, and 
from positions off the line of sight. We show that at time t = r/2^ 2 {r)c the fireball appears on the 
sky as a narrow ring of radius r/ , j(r) (This conclusion is independent of the details of electron and 
magnetic field distribution). Our conclusions and their implications are summarized in §3. 



2. Emission from a relativistic expanding blast wave 

Let us consider a strong, spherical, ultra-relativistic shock wave expanding into an ambient 
medium of uniform density. For a shock propagating with a Lorentz factor r> 1, the Lorentz 
factor, number density and energy density of shocked fluid at the shock discontinuity are 
7 = r/v2, n = 4qri{ and e = 4j 2 nim p c 2 respectively, where rij is the number density ahead of 
the shock (e.g. |Blandford &; McKee 1976| ). Observations indicate that the fraction of blast wave 
energy carried by electrons is not large, ~ 0.1, and that the electron cooling time is long compared 
to the dynamical time, i.e. to the blast wave expansion time ( [Waxman 1997a ,b). This implies that 



the energy lost to radiation is small, and that the blast wave energy is approximately constant 
(so called "adiabatic blast wave"). Since the total energy in shocked particles is proportional to 
njr 3 7(e/n), conservation of energy implies 

/ r \ ~ 3 / 2 

7 = 7o — , (1) 



where ro is some fiducial radius. 

In fireball models of GRB afterglow it is assumed that the fractions of energy carried by 
magnetic field and by electrons are time independent. Under this assumption the magnetic field 
B and the characteristic electron Lorentz factor 7 e (in the shell frame) scale as B cx 7 e oc 7. The 
observed radiation is produced by synchrotron emission of the shock accelerated electrons. The 
characteristic synchrotron frequency in the shell frame scales as v oc j 2 B oc 7 3 . The number of 
photons emitted as the shock propagates a distance dr may be obtained as follows. The number 
of radiating electrons scales as r 3 , and the number of photons each electron emits per unit time 
in the shell frame is proportional to the magnetic field B. The time in the shell frame over which 
the shock propagates a distance dr is dr/^c. Thus, the number of photons emitted scales as 
dN/dr oc r 5 B/j. Using (|l|) we therefore have 

„ = V0 (L)-"\ WiV, (2) 



r J dr \r 
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where L is a constant with dimensions of length. 

The shock heated gas expands relativistically in its rest frame. Since the time, measured in 
the rest frame, for the shock to expand to radius r is ~ r/jc, relativistic expansion in the shell 
frame implies that the rest frame thickness of the shell of shock heated gas is ~ r/7. Thus, in the 
observer frame most of the shocked gas, and most of the blast wave energy, are concentrated in 
a shell of thickness Ar = C, r /l 2 i where £ is some constant. We can obtain an estimate of Q by 
assuming that the density in the shocked shell is uniform. In this case conservation of particle 



number implies Airr 2 'yAr4'yn = 4-/rr 3 n/3, i.e. £ = 1/12. In the self-similar solutions of ( Blandford 



fc McKee 1976| ), which give the spatial dependence of the hydrodynamic variables, 90% of the 



energy is concentrated in a shell of thickness corresponding to £ = 1/7. It is implicitly assumed in 
fireball afterglow models that the fractions of energy carried by electrons and magnetic field vary 
behind the shock over a length scale comparable to the scale for changes in the hydrodynamic 
variables, Ar. This is a reasonable assumption, since the synchrotron cooling time of the electrons 
is longer than the dynamical time, i.e. the characteristic expansion time. If the energy fractions 
vary behind the shock on a scale much shorter than Ar, due to some non hydrodynamic process 
which is not accounted for, a new length scale would be introduced into the problem and the 
scalings @, on which the fireball afterglow model relies, will no longer hold. Since the details 
of the spatial dependence of the electron and magnetic field energy fractions are not known, we 
will assume below that radiation is produced within a homogeneous shell of width Ar behind 
the shock, and will derive results for different values of £. For clarity we first assume that at a 
fixed time all emitted photons have the same energy (in the fireball frame). We then discuss 
modifications expected due to a finite frequency range of emitted photons. 

In the shell frame photons are emitted isotropically. Consider, therefore, the appearance 
of an isotropic distribution of photons of frequency v in a frame moving with Lorentz factor 7 
with respect to the frame where the distribution is isotropic. We refer to the frame where the 
distribution is isotropic as the "rest frame" . Denoting with primes quantities measured in the 
moving frame, the rest frame and moving frame frequencies are related by 

v' = 7(1 + /3coe0)i/, (3) 

where (3 = (1 — I/7 2 ) 1 / 2 and 9 is the (rest frame) angle between the photon momentum and the 
direction of motion of the moving frame. The angle measured in the moving frame is 

„; sin 9 , AS 

7(p + cos 9) 

The fraction df of photons with frequencies in the range v' to v' + dv 1 is 

df 1 . a d9 1 . , , , ,_. 
^ = 2 Sm ^7 = 2^ ^„<-<W (5) 

Here > U Ln = 7(1 - P)V, z4 ax = 7(1 + P)V. 



- 5 - 



Let us now consider photons observed by a distant observer in a frequency range v' to v' + dv' . 
The number of photons produced by the fireball shell at radius r, with observed frequency in the 
range v' to u' + dv' , is given by (Q) and (||) 



d 2 N 
dv'dr 



L 



-i 



2^{r)v(r)rl 2L^ Q v \r 



(6) 



where we have approximated j3 = 1 in (||), as we are interested in the limit 7 > 1. Any photon 
which arrives at a distant observer must be emitted from the fireball in a direction parallel to the 
source-observer line of sight. A photon emitted at radius r, with frequency v{r) in the shell frame, 
is observed with frequency v 1 provided, therefore, it is emitted from the fireball on a line emerging 
from the explosion center and making an angle 9', given by @ and (|3|), with the source-observer 
line of sight. Such a photon is delayed with respect to photons emitted on the line of sight by 



At e = (1 - cos 1 



P 



El 

c 



7(1 + 0): 



1 



2 7 2 c 



27- 



1 



(7) 



The total delay, i.e. the delay with respect to photons emitted from the center of the explosion 
r = 0, is given by the sum of Ate and At r = t e — r/c, the difference between light travel time to 
radius r and shock expansion time to radius r, t e = f£ (1 — l/r 2 ) _1//2 dr/c. For the scaling (jl|) we 
have, for 7>1 and using T = 2 1 / 2 7, At r = r/167 2 c. Thus, the delay of photons emitted from the 
shock front at radius r, and observed with frequency v 1 ', with respect to photons emitted from the 
center of the explosion r = 0, is 



t{u', r) = At e (u', r) + r/l^ 2 c. 



(8) 



Since photons are emitted uniformly from a shell of finite thickness Qr/^ 2 c behind the shock, 
the arrival times of photons emitted at radius r and observed with frequency v' are uniformly 
distributed (for 7»1) over the range t = t(v', r) to t = t(z/, r) + 7£ , where = (r/~/ 2 c. 

Using (|6|-||) we have numerically calculated the intensity (energy flux per unit frequency) 
f v = v'd 2 N/dv'dt as a function of time. Results are shown in Fig. 1 for different values of the 
shell thickness, Q = 1/4, 1/16, and 1/64. Fig. 1 presents the time dependent flux at v 1 = 70^0- 
The intensity at other frequencies is obtained using the scaling 



, , v. 



-2/3 ' 

- ! t 



(9) 



For the shell width expected from fireball hydrodynamics, Q ~ 1/10, the intensity peaks at 
frequency v' = j{r)u(r) at a time t — r /A^ 2 {r)c. The numeric coefficient in this relation, 1/4, 
is somewhat smaller than the commonly used value, 1/2. However, as explained above, since 
the spatial dependence of the fractions of energy carried by electrons and magnetic field is not 
known, the effective shell width Q can not be determined exactly (as it is not determined by the 
hydrodynamics only), and therefore the value of the coefficient in the relation t ~ r/47 2 (r)c can 
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not be determined exactly. It is nevertheless clear that the intensity peaks at v 1 = ^(r)v(r) at a 
delay significantly larger than t = r/167 2 (r)c. 

The results presented in Fig. 1 were obtained under the assumption that at a fixed shell 
radius all emitted photons have the same frequency v(r) (in the fireball frame). In fireball 
afterglow models the emission peaks at v(r), and extends to high frequency v 3> v(r) due to 
power-law energy distribution of electrons, and to low frequency v <C u(r) due to the low energy 
tail of synchrotron emission, f(v)/f[v(r)] ~ [u /^(r)] 1//3 . These tails dominate the intensity at 
observed frequency v' = j{r)u(r) at delays much smaller and much larger than the time at which 
the intensity peaks at u', t <C r/^/ 2 c and t 3> r/7 2 c. However, they are not important near the 
peak, t ~ r/7 2 c. The details of the frequency dependence of the emission at v ~ v(r) will affect 
the behavior at t ~ r/j 2 c. Let us assume that photons are emitted at the shell frame at two 
different frequencies, v(r) as given in (Q) and u(r) = xu(r). It is straight forward to show that the 
intensity fv(v',t) due to emission at v{r) is related to that of the intensity /„(i/,t) due to emission 
at u(r) by f u (i/',t) oc f v {v',x~ 2 ^t). Thus, if most of the energy is emitted in the shell frame over 
a frequency range 8v(r)/v{r) ~ 1, the frequency spread would introduce a spread in arrival time 
of photons of frequency v\ 5t/t p ~ 1 where t p is the time at which the intensity peaks under the 
assumption that all photons are emitted at the same frequency v{r). From Fig. 1, the spread in 
arrival time due to the relativistic expansion is 5t/t p ~ 2. Thus, the effect of emission over a finite 
frequency range would not be significant, provided most of the energy is emitted over a frequency 
range 8v{r)/v(r) <^ 1. 

The flux observed at a fixed time originates from points at a range of distances transverse to 
the line of sight. Using the equations derived above we have numerically calculated the fraction 
of flux contributed from rings of radii h— h + dh around the line of sight, as function of h. The 
results are shown in Fig. 2 for several £ values. At time t = r/2^/ 2 (r)c the flux originates from a 
ring of outer radius h ~ r/j(r) and width Ah/h ~ 0.1. Note that this result is independent of the 
spectral distribution of the emission in the shell frame. The appearance of the fireball as a narrow 
ring can be qualitatively understood from Fig. 3. Radiation produced at ro and observed at 
*o = r o/^1 2 ( r o) c originates from an angle ~ l/7o> corresponding to h = ro/70 Contribution to the 
emission at smaller h is due to emission from r < ro and r > tq. However, the contribution from 
smaller radii is suppressed due to the fact that at smaller radii the radiation seen at a given time 
originates from larger angles, and the relativistic beaming suppresses the emission from angles 
> I/7. The contribution from larger radii is small since the emissivity decreases as the fireball 
decelerates. 



3. Conclusions 

We have presented a detailed analysis of the implications of the relativistic hydrodynamics 
of expanding blast waves to GRB afterglow observations. We have shown that the afterglow 
intensity peaks at observed frequency v 1 = j(r)v(r), where v(r) is the frequency at which the 
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emission peaks in the fireball frame at radius r, at time t ~ r/47 2 (r)c. The exact value of the 
numerical coefficient depends on the effective thickness of the radiating shell, C r /7 2 - The value 1/4 
is obtained for £ ~ 1/10, the value implied by fireball hydrodynamics (see Fig. 1). The numeric 
coefficient, 1/4, is somewhat smaller than the commonly used value, 1/2. However, it should be 
kept in mind that since the spatial dependence of the fractions of energy carried by electrons and 
magnetic field is not known, the effective shell width £ can not be determined exactly, and the 
numerical coefficient may be somewhat larger or smaller then 1/4. The amplitude of the peak 
intensity agrees with the estimate of QWaxman 1997a| ,b). For the intensity normalization chosen 



in Fig. 1, the peak intensity estimated following ( Waxman 1997a ,b) is 2/3, close to the numerical 



results obtained here. This is due to the fact that most photons of frequency v' = j(r)v(r) arrive 



over a time scale t ~ r/2j (r)c, as assumed in (Waxman 1997a, b) 



At time t the fireball appears on the sky as a narrow ring of radius h = r/j(r) and width 
Ah/h ~ 0.1, where r and t are related by t = r/2j 2 (r)c (see Fig. 3). The apparent size 
h = r/"i(r) is smaller by a factor of (16/2) 5 / 8 ~ 4 compared to that obtained by using the relation 
t = r/167 2 (r)c (this follows from |l]]). The smaller size implies that diffractive scintillation is likely 
to modulate the afterglow radio flux ( |Goodman 1997 ), and that significant modification of the 



optical light curve due to microlensing is possible on day time scale ( Loeb fc Perna 1997 ). The 
narrowness of the emission ring would affect the predictions for both microlensing and scintillation 
(This has been taken into account in the lensing calculations of [ [Loeb fc Perna 1997 ] . More 
detailed calculations are required for modulation by scintillation). 

The results presented here are valid for highly relativistic fireballs, 7>1. For Lorentz factors 
7 — 1 ~ 1, the emission ring would be wider. Finally, we note that although we have implicitly 
assumed spherical symmetry throughout the paper, our results are valid for a fireball which is a 
cone of finite opening angle 9 as long as 7 > 1/9. This is due to the fact that most of the observed 
emission originates from a cone of opening angle I/7. 
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Fig. 1. — The observed intensity (energy flux per unit frequency) at observed frequency 
v 1 = j(ro)v(ro) as a function of time for various fireball shell widths. The width expected from 
hydrodynamic considerations corresponds to £ ~ 1/10. v(tq) is the frequency with which photons 
are emitted in the shell frame at r = vq. The normalization of the intensity is determined by 
measuring time in units of ro/7 2 (ro)c, and by setting the dimensional coefficient in (Q) to L = 1. 
The vertical solid line denotes the time t = ro/167 2 (ro)c. 
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Fig. 2. — The fraction of flux observed at time t = ro/2'y 2 (ro)c which originates from different 
distances h transverse to the line of sight, for various fireball shell widths Qr/^ 2 {r). 
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Fig. 3. — The fraction of flux obserevd at time t = ro/27 2 (ro)c and produced over a radii range 
dr, df /dr, plotted as a function of r with h(r), the distance transverse to the line of sight from 
which radiation emitted at r arrives at the observer at time t, for a shell of thickness £ — > 0. 



